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Motivated by recent experimental work on FeCrAs, we study the magnetic properties of classical 
spin models on the hexagonal Fe2P/ZrNiAl structure (space group P&2m). When both transition 
metal/rare earth sites carry magnetic moments, one has alternating distorted kagome and triangular 
lattice layers. Each point of the triangular lattice lies at the centre of a distorted hexagon and 
actually corresponds to three isolated magnetic atoms: a "trimer". We show that a simple model 
consisting of antiferromagnetically correlated Heisenberg spins on the frustrated distorted kagome 
lattice coupling either ferromagnetically or antiferromagnetically to Heisenberg spins at the average 
position of a trimer, leads to a rich phase diagram, including the possible selection of a (i, i,0) 
magnetic order as may have been seen in FeCrAs. 

PACS numbers: 75.10.Hk,75.40.Mg,64.60.Ej 



INTRODUCTION 

Materials featuring magnetic moments situated on 
corner-shared triangles or tetrahedra provide an exciting 
domain for the discovery of exotic physics. Recently one 
class of such geometrically frustrated magnets (GFM), 
spin-ice materials, garnered attention due to the possibil- 
ity of interpreting the excitations from the ground states 
as effectively deconfined magnetic monopoles [l[ . Several 
experimental groups have claimed Q to have observed 
these so-called monopoles, and have pointed out that 
they would likely pass several of the experimental tests 
for the existence of magnetic monopoles. Of particular 
interest was that these fractional spin excitations appear 
to be deconfined or essentially free in three dimensions. [l| 
Recently we have shown that spin ice physics should also 
arise on the three-dimensional trillium and hypcr- 
kagome 0] lattices. 

In a second class of GFM, Heisenberg spin systems, a 
long sought goal has been the realization of a truly three- 
dimensional quantum spin liquid. Indeed, the belief in 
the existence of such a state strongly influenced the orig- 
inal resonating valence bond model of high temperature 
superconductivity Therefore it was with great excite- 
ment that the community welcomed the Mott insulating 
Na4lr308 which features antiferromagnetically correlated 
spin 1/2 magnetic moments and does not magnetically 
order to very low temperatures Q. These phenomena are 
well captured by a simple classical Monte Carlo treat- 
ment of the Heisenberg model on the three dimensional 
corner-shared equilateral triangle hyperkagome lattice^], 
which additionally predicted the onset of a finite temper- 
ature nematic transition. Early quantum treatments of 
the physics of this material have argued that at least 
over a finite temperature range one has indeed found 
a quantum spin hquidli-[ii| or at least a valence bond 



crystal. [12 1 

It has been considerably more difficult to experimen- 
tally realize such exotic physics on the simpler two- 
dimensional kagome net, perhaps because materials are 
inherently three-dimensional. Indeed, after decades of 
searching, the first experimental realization of antiferro- 
magnctic correlations on a structurally perfect spin 1/2 
kagome lattice, ZnCu3(OH)6Cl2, 13|,ll4| appeared only in 
2005, followed soon after by the related ZnCu3(OD)KCl2 
which shows evidence of a valence bond solid state [l5| . 
However, even in these materials questions remain about 
the realization of a purely kagome lattice, given the re- 
cently measured non-magnetic site disorder of the order 
of 6% within the magnetic Cu plane 0. 

Interestingly, a large class of magnetic materials in- 
cluding some (non-superconducting) iron pnictides and 
the langasites form in a distorted kagome structure 
with a nearest neighbour topology equivalent to that 
of the kagome lattice. Among these, Nd langasite 
(Nd3Ga5SiOi4]has previously been studied as a spin liq- 
uid candidate [17|, while FeCrAs shows intriguing ther- 
modynamic and transport signatures. [31 Magnetically, 
FeCrAs has been seen[l9| to exhibit a -s/S x -s/S magnetic 
ordering with wavevector of (i, i, 0) below 125K. 

While the nature of the spin-^ quantum Heisenberg 
model on the kagome lattice remains an important open 
problem j20j. the classical antiferromagnetic Heisenberg 
model has been shown [2l] - [25l | to have degenerate so- 
called (7 = and -y/S x -y/S ground states, and exhibit an 
order by disorder transition to select the latter at T — )■ 0, 
resulting in a local nematic orderin g w ith a diverging cor- 
relation length in the limit T — 0[21|. Why then would 
FeCrAs choose to magnetically order at 125 K into a 
(i,i,0) ground state? 

Could it be that nearest and next nearest neighbor ex- 
change plus Dzyaloshinskyy-Moriya interaction terms are 
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all necessary to a description of the physics on the dis- 
torted kagome lattice as has been studied by Gondek et 
aZ. [2^ for XY spins, which gives rise to interesting mag- 
netic structures, including one with a (^,5) propagation 
vector? The nearly isotropic resistivity and magnetic sus- 
ceptibility, and the existence of the phase transition at 
finite temperatures speak to the possibility that we have 
a truly three dimensional material, naturally leading us 
to ask what the effect of a magnetic coupling between 
frustrated distorted kagome planes would lead to. 

Curiously, the addition of an interlayer coupling be- 
tween stacked kagome planes to an antiferromagnetic 
Heisenberg model, and its role in leading to possible mag- 
netic order, appears not to have received much atten- 
tion. The spin-^ quantum antiferromagnetic Heisenberg 
model has been treated on the stacked kagome lattice 
by D. SchmalfuB et aZ.|27| using a spin-rotation-invariant 
Green's function method (supported by classical spin- 
wave calculations), finding the system to remain "short- 
range ordered independent of the strength and sign of the 
interlayer coupling." For classical spins, M. Zelli et al. 28 1 
used a short-time dynamics Monte Carlo method, scal- 
ing analysis and Binder cumulants to study the ordering 
properties of a generalized antiferromagnetic Heisenberg 
model on a stacked triangular lattice. A limiting case 
of this model is the classical Heisenberg model on the 
stacked kagome lattice for which they found no evidence 
of long range order. 

However, when the couplings between kagome planes 
are mediated by other magnetic sites, a dramatically dif- 
ferent result can be expected. FeCrAs is an example of 
such coupled kagome systems in which both Fe and Cr are 
magnetic ions, and the dominant magnetic interaction is 
believed to be an antiferromagnetic Heisenberg coupling 
between local moments on the Cr sites. The Cr sites form 
stacked distorted kagome planes. Three nearest Fe ions 
form a trimer (triangle) , which can be viewed as an effec- 
tively single magnetic site. These effective Fe molecules 
then form triangular planes between Cr kagome planes. 
Therefore. FcCrAs can be viewed as alternating distorted 
kagome and triangular planes as shown in Fig. [TJ A mi- 
croscopic derivation of this physics, and the validity of 
our approximation will be presented elsewhere. [2^ 

In this paper, we propose a minimal microscopic model 
that captures the magnetic ordering in FeCrAs at finite 
temperatures. Using MC simulation, we find that any 
small finite coupling between the Fe and Cr sites leads 
to a coplanar a/S x ^/S magnetic order consistent with 
existing neutron scattering data on FeCrAs. As the 
coupling between Fe and Cr increases, there appear two 
finite temperature transitions. At high temperatures one 
has a transition to a ferri(ferro)magnetic state followed 
at low temperatures by a canting of each spin to realize a 
component of each spin which adopts a y/S x order of 
some fraction of the spin. At larger interplane couplings 
still this gives way to ferri(ferro)magnetic order, as is 



commonly seen in other members of this family. This 
implies that we may have a large family of previously 
unrecognized frustrated magnets, and should have im- 
plications for uniaxial pressure experiments on FeCrAs, 
which can be tested by a future experiments. 



GENERAL PROPERTIES AND LATTICE 
STRUCTURE 

Lattice structure 

Materials forming in the hexagonal Fe2P/ZrNiAl struc- 
ture with space group P62ml30| feature two site types 
potentially occupied by magnetic atoms labeled 3/ and 
3g as seen in Table U and shown in Fig. [TJ When 
in a given layer the value of i < < | (i = 1 or 
2), the magnetic atoms belonging to that layer are ar- 
ranged in a distorted kagome net, topologically equiva- 
lent to the kagome net for nearest neighbour couplings. 
Thus antiferromagnetic interactions or ferromagnetic in- 
teractions with a strong single ion anisotropy would be 
strongly frustrated, although away from the special point 
Xi = ^, which corresponds to the kagome lattice, the 
lattice lacks inversion symmetry so that Dzyaloshinkii- 
Moriya terms can arise. Alternately, when < Xi < ^ or 
I < < 1, nearest neighbour bonds between magnetic 
sites form isolated triangles about the centers of the dis- 
torted hexagons in the layers directly above and below. 
The special limit Xi = | in a layer corresponds to the tri- 
angular lattice. For simplicity in our model below we will 
take the limit xi — which treats the isolated trimers as 
all being at a single point (0,0,0) of the hexagonal lattice. 



Site type 


Atomic position 


3/ 


("i",-^\"i",0), ("i",^\"i",0), {-xia,0,0) 


3g 


/ X2a V3x2a c\ / X2a ^Zx^a c\ ( n. c\ 


Name 


Lattice vector 




/ a v^a r\\ 




(f,4^,o) 




(0,0,c) 



TABLE L Atomic positions of potentially magnetic atoms 
{3f,3g) in the Fe2P/ZrNiAl structure, and hexagonal lattice 
vectors {Ai, A2, A3}. Generically xi < ^ which leads a layer 
of isolated "trimers", and | < 2:2 < | which leads to a dis- 
torted kagome lattice of corner-shared equilateral triangles 
except at the special point X2 ~ ^ which corresponds to a 
perfect kagome lattice. 



Relevant magnetic materials 

A large number of magnetic materials form in 
the hexagonal Fc2P/ZrNiAl structure. It has been 
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estimated [26| that almost 30% of all rare-earth- 
transition metal compounds RTX crystallize in this 
form, with the rare earth metal generally sitting in 
a distorted kagome site. Among these, antiferromag- 
netic order is seen in RAgSi for R=|Dy,Ho,Er,Tb j |31| 
and RAgGe for R={Gd,Er,Tb,Dy,Ho}0|, which are 
metallic materials for which the RKKY interaction 
and strong crystalline electric fields have been thought 
to be relevant. Antiferromagnetic Curie-Weiss offsets 
have likewise been seen in TbPtIn and TmAgGe, which 
have been argued to be describable in terms of a triple 
coplanar Ising model with large moments consistent 
with R^+.js^ A large number of TT'X materials also 
form in this structure, some showing evidence of ferro- 
(or possibly ferri-) magnetic order, others showing 
evidence of antiferromagnetic order. With T=Ni and 
X=P, T'={Mo,W,Fe,Co,Cr[33} are known as are 
NiMnAsllil and MnTiP. Known arsenides feature TT' 
={CrPd,FeV,MnRu,CoCr}[35|. KnownH^ silicides 
include TT'={TiMn,ZrRu,NbCr,NbMn}, and knownjsil 
germanides include TT'={NbCr,HgFe,NbMn,TiCo}. 
Amongst the langasites one has at least two magnetic 
mcmbersji^, Nd3Ga5SiOi4 and Pr3Ga5SiOi4, although 
in these materials the trimer lattice is occupied by Ga 
which is unlikely to have a magnetic moment. 



the magnetic interaction can be described as an antifer- 
romagnetic Heisenberg model between local moments on 
the Cr sites (/i ~ 1.22^b[i^) that form distorted kagome 
planes as shown in Fig. [IJa). In between the Cr planes 
lie Fe atom trimers which have a small moment (if any 
fi < (0.1 ± 0.03)/z_b[i^) and lie equidistant directly be- 
low and above the centre of the distorted hexagons of 
Cr in adjacent planes. Since according to the Mermin 
Wagner theorem one does not expect lon g ra nge order to 
arise in a purely two dimensional system [37|, we believe 
that at finite temperature an interplay between frustra- 
tion and dimensionality may give rise to this behavior. 
Namely, we posit a nearest neighbour coupling between 
the distorted kagome Cr planes and the Fe trimers that 
lie between. A competing candidate for the origin of the 
dimensional crossover to allow a finite temperature phase 
transition if the Fe sites do not hold magnetic moments 
would be a dipolar interaction between neighboring Cr 
planes, although such a picture would need to account 
for the surprisingly high ordering temperature, given that 
dipolar interactions are estimated to be only of the order 
1.27 X IQ-^K here. 



FeCrAs 

In FeCrAs it is believed[l^ that there is perfect or- 
der of the Cr atoms in the 3^ sites, and the Fe atoms 
in the 3/ sites. Swainson et al.^^ claim that for this 
material xi = 0.760(2), X2 = 0.436(2), a = 6.06761, and 
c = 3.6570A, although these detailed values do not enter 
into our calculations of the heat capacity and magnetic 
susceptibility of the model presented below (which have 
additionally taken the xi — >■ (or equivalently 1) limit). 
For these values, we note that the nearest neighbour Cr- 
Cr bonds form a distorted kagome lattice, while the Fe-Fe 
bond distance is quite short. 

FeCrAs is an unusual magnetically frustrated iron 
pnictide that has recently garnered interest as a possi- 
ble novel non- Fermi liquid [18[. FeCrAs shows 18j Fermi 
liquid behaviour with a specific heat coefficient 7 of 
30mJ/molK^ which indicates that it may be metallic 
but this contrasts with a resistivity which increases al- 
most isotropically from &00K to 80miir[l^. Apart from 
a feature near a magnetic ordering transition in the a-axis 
resistivity, this increase at low temperatures is mono- 
tonic and below lOK appears to follow a weak power 
law Sp ~ — (q; ^ 0.6 — 0.7) along both the a and 
c-axesfisj rather than the hiT one might expect from 
Kondo physics. 

At low temperature (T ^ 125K) FeCrAs shows a mag- 
netic transition with an ordering wave vector of fc = 
(5, 5; 0)[l3, E^. From the susceptibility measurement. 



a) 



b) 




FIG. 1: (Color online) a) The FeCrAs lattice with (blue) cir- 
cles at Cr (3g) sites joined to nearest neighbors and (red) 
circles at Fe (3f) sites, b) The groups of three Fe sites are 
taken to act together as a single Fe trimer. 



4 



Binder 



,0.5- 



• • 

V3xV3 Canted" 
W3 X V3 



♦ 

Ferrimagnetic 



0.5 1 1.5 

yjj (Withjj + ij2i = 1) 



2.5 



FIG. 2: (Color online) The phase diagram of FeCrAs as de- 
termined from peaks of heat capacity vs. temperature plots 
of Monte Carlo results at L = 6. Also included are the transi- 
tion temperatures obtained from the Binder cumulant. These 
temperatures lie only slightly above the results obtained from 
the heat capacity data, indicating that the ferrimagnetic tran- 
sition is not an artifact of the finite lattice size. The dotted 
line is a guide for the eyes. For J2 < the ferrimagnetic 
region is simply replaced by a ferromagnetic region. 



MODEL 

We consider a classical antiferromagnetic Heisenberg 
model between nearest neighbour Cr spins as has been 
previously studied in the context of the Nd langasites 
by Robert et al.^^. For simplicity, to address the in- 
teresting dimensional crossover question, we assume that 
the iron atom trimers act as a single classical spin (see 
Fig. [Ijb)) and that there is a Heisenberg coupling be- 
tween each of these Fe trimers and Cr sites. The model 
is written: 



H 



■ Cr 



ilk) 



' Cr 



• Fe 



(1) 



where Ji > is an antiferromagnetic coupling and J2 
may be antiferromagnetic (J2 > 0) or ferromagnetic 
(J2 < 0). For this simple model, wc assume that the 
spin coupling between the Cr and Fc sites is the same for 
the two inequivalcnt (but approximately equal) neigh- 
bor distances from the trimers. Thus each Cr site has 4 
Fe neighbors with the same magnetic interaction in this 
model. Noting that the magnitude of any non-zero net 
spin on the Fe trimer can be (at least energetically) ab- 
sorbed into the definition of J2, we proceed to describe 
the general phase diagram as ^ is varied from to inf 
keeping Ji + IJ2I ~ 1 and unit spins on each site. 



Ground State Energy 

The ground state energy for various values of Ji and 
J2 was found using numerical minimization. As we in- 
crease J2 from zero, we see that initially (i.e. -j^ = 0) wc 
recover the known ground states of the Heisenberg model 



on the kagome lattice, which feature 120° rotated spins 
on each of the corners of the equilateral triangles of the 
Cr sites, and a randomly oriented spin at the Fe site. At 
the opposite extreme, when ^ — >■ oo(— 00), we have a 
ferrimagnetic(ferromagnetic) structure with Cr spins in 
the plane aligned antiparallel(parallel) to Fe spins, but 
parallel to each other. In between the spin structure has 
components of each structure, at each Cr site adopting a 
form, 



S = (cos((?!))si2o,sin((/))). 



(2) 



That is, for any three Cr spins forming a trian- 
gle, one can choose (j) such that the spins are given 
by (cos(,/>),O,sin(0)), (-£2^, , sin(0)), and 

/ cos{4>) \/3 cos(</)) 



,sin(0)). 



The 



2 

Fc 



spins are 



then 



V 2 ' 2 

opposite(parallel) to the sums of their nearest two spins. 
Since the Fe spins lie above and below the centers of 
the distorted hexagons of the lattice, this implies that 
each hexagon has only two spin types around it. That 
is, the so-called -\/3 x -\/3 ordering of the kagome lattice 
is selected over the q = state by the presence of 
the Fe atoms. This explains why there is a (■|,-i,0) 
periodicity in the magnetic lattice as reported by Wu et 
al. [l8| , as it takes three unit cells in each of the in-plane 
lattice vector directions for the spin structure to repeat. 
Next, let us compute the energy of this ground state. 
For completeness, the directions of the Fc atoms can 
as -«gn(J2) ^cos^^ ^7(^\ 2sin(0)), 



s(0) 



V3 cos(0) 



2sin((/))), 



and 



be taken 

-sgn(J2) 
^l+3sin2(0) 

^^ii^(-cos(0),O,2sin(</))). Carrying out the 

dot products wc find that the energy per spin (noting 
that there are 3 Cr and 1 Fe sites per unit cell in our 
model) is then given by, 



3( Ji(^sin2(0) 



)-|J2K/l + 3sin^(0) 



(3) 



Minimization of e with respect to (j) yields possible solu- 
tions of sin(0) = 0, cos(0) 



Oand 1^ = y'l + 3sin2(,/)) 

(the latter being valid only for 1 < -^^^ < 2). The former 
(0 = 0), corresponds to the spins maintaining a purely 
120° rotated spin structure in the Cr plane, with the 
spins on the Fe sites coplanar, and opposite(parallel) to 
the spin sum of their nearest Cr sites. The energy of 
this state can be expressed as — — 3|J2|- This nicely 

reproduces our minimized energies until = 1. At 
this point, the Cr spins begin to cant upwards uniformly 
from the 120° rotated spin structure, while the Fe spins 
begin to tilt downward (upward if ferromagnetic). The 

3j2 

energy of this spin structure is given by — 3Ji — -jj^ , 
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which minimizes the energy until ■'-j^ = 2, at which point 
(j) = ^. Beyond this point, the spin structure is ferri- 
magnetic (ferromagnetic), with an energy per unit spin of 
3Ji — 6| J2I which corresponds to all Cr spins parallel, and 
Fe spins antiparallel(parallel). 



indicating that the Binder results are consistent with the 
thermodynamic limit of our results. 



HEAT CAPACITY 



MONTE CARLO SIMULATION 

Monte Carlo simulations were carried out using the 
Metropolis algorithm with periodic boundary conditions 
on a lattice with hexagonal unit cells with a total side 
length of L unit cells containing AL^ sites in total (| be- 
ing Fe trimers and | being Cr atoms) . The bulk of the 
analysis was done using a lattice size of L — 6 which cor- 
responds to 864 spin sites. At each temperature 2 x 10^ 
Monte Carlo steps were used to equilibrate the system 
and a further 2 x 10^ were used to calculate the averages 
of physical quantities, where one Monte Carlo step was 
taken to, on average, attempt one update per site. Error 
bars as reported correspond to the standard deviation of 
our averages over four independent trials. Each update 
attempt has a randomly chosen magnitude, 6, and direc- 
tion relative to the initial spin direction. The value of 
S was chosen so that around 50% of the attempted spin 
updates are accepted for any combination of Ji and J2 
at all temperatures nil- 

In order to take into account any finite size scaling 
effects the Binder cumulant [sij was calculated for several 
values of 4^. The Binder cumulant is defined as 



where m is the magnetization of the lattice. For the fer- 
rimagnetic transition we have defined m = ^ 

SfeLi s^'^ which clearly demonstrates that in the region 
featuring two transitions, the high temperature transi- 
tion indicates the onset of ferrimagnetic order. The cu- 
mulant has a critical temperature that is independent of 
the size of the lattice and is defined by the intersection 
point of the data acquired from simulations using differ- 
ent lattice sizes. The critical temperature obtained from 
the Binder cumulant is important as it defines a size in- 
variant transition temperature for the lattice. This allows 
us to compare the transition temperatures extracted from 
the L = 6 heat capacities with results that should hold to 
the thermodynamic limit. Fig. [2] shows that the transi- 
tion obtained from the Binder cumulant differs from the 
values obtained from the heat capacity measurement by 
only a small positive amount. This small difference in- 
dicates that the transition temperatures obtained from 
the heat capacity were not caused by finite size effects 
and are reasonably close to the values that would be ob- 
tained in the thermodynamic limit. The L ~ 9 (not 
shown) results are seen to lie closer to the Binder results. 



To determine the magnetic phase diagram of this 
model we have calculated the heat capacity func- 
tion of temperature as a function of Ji and J2. From 
peaks in the heat capacity, we can extract the order- 
ing temperature for various Ji (with Ji -I- IJ2I — 1). 
This allows us to write an approximate phase diagram 
delineated by the transition temperatures involved, as 
shown in Fig. [2] In the limit where > 2 one has 
a ferri(ferro)magnetic form with spins in the Cr plane 
aligned and antialigned(aligned) to the spins in the Fe 
trimer. When 1 < 1^ < 2 (Fig. [SKa)) the Cr spins ac- 
quire both a component which has 120° rotated spins and 
a component which tilts upward relative to the Cr plane, 
so the spin structure is canted. One appears to have two 
clear transitions. In the last region where < < 1 
(Fig. [31[b)), all eases have a 120° rotated spin structure 
in the Cr plane, and coplanar spins on the Fe trimers 
whose sum is opposite(parallel) to the spin sum of the 
nearest Cr sites. 

Representative plots of the heat capacity in units of 
the Boltzmann constant per atom as a function of tem- 
perature are shown in Fig. [3l As these are classical 
Heisenberg spins subject to the constraint of unit magni- 
tude, one expects from the equipartition of energy theo- 
rem that at low temperatures the heat capacity will have 
two quadratic degrees of freedom per site, leading to ^ks 
each, for a total Cy = fc^. At low temperature in real 
quantum materials this cannot be physical, however, it is 
hoped that the higher temperature classical phase tran- 
sitions may remain relevant. In this sense we see that 
as we turn on the coupling between the Cr atoms (from 
Jl = 0), the transition temperature decreases monotoni- 
cally until the point = 2 below which two transitions 
(see Fig. EJa)) occur as a function of temperature. This 
coincides with the transition from a ferrimagnetic ground 
state to establishing a ground state with a component 
which rotates on each triangle of the lattice (0 7^ 0). At 
the high temperature transition, the ferrimagnetic orien- 
tation of the spins sets in, with the coplanar 120° com- 
ponent attaining order at still lower temperatures. From 
the low temperature transition (see Fig. ^h)) the tran- 
sition temperature rises a little for the next two values 
of ^ = 1 and -j^ = I before decreasing monotonically 
towards zero temperature for the uncoupled layer case. 

The heat capacity of the kagome case (^^ = 0, Fig. 
[3ljb)) does not approach 1 but instead approaches | as 
the temperature approaches zero as the curves for other 
values of Ji do because of its normalization. Energeti- 
cally it is as if the Fe trimers have zero spin as they no 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Temperature [Jj + J^] 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Temperature [J^ + J^] 

FIG. 3: (Color online) The heat capacity of the system 
L = 6 for: a) 1 < ^ < 2, and b) < 1. The values of 
are labeled on each plot. 



at 
£2 

Ji 



longer have any contribution to the heat capacity of the 
lattice, yet we have included them in the normalization. 
No ordering transition is seen, although a weak peak is 
observed in the heat capacity at low temperature consis- 
tent with previous work done by Zhitomirsky et aL[40| 
and others 21 1. 



MAGNETIC SUSCEPTIBILITY 

To guide the development of a realistic model of 
FeCrAs, we have performed calculations of the static 
magnetic susceptibility as a function of the temperature, 
^(T) = jpj, ^ < sl{T) ■ s'j{T) >. It is anticipated that 

id 

such studies of classical magnetism on a frustrated lat- 
tice may show qualitatively similar results to a quantum 
version of our model, allowing us to estimate realistic pa- 
rameters for Ji and J2. The bulk of our analysis to date 
has been carried out with J2 > 0. The magnetic suscep- 
tibility for 7^ < has been found to follow similar trends 
to that of the ^ > case with the main difference being 




0.5 1 1.5 2 2.5 
Temperature [Jj+J2] 



FIG. 4: (Color online) The inverse magnetic susceptibility 
as a function of temperature. For ^ = 9.0 to 1.5, the in- 
verse susceptibility quickly falls towards 0. For ^ — 1 this 
quick descent is followed by a brief low temperature upturn. 
Qualitatively different physics is seen for -jj- < |, where one 
sees a high temperature linear dependence, followed by a low 
temperature downturn (likely due to the paramagnetic contri- 
butions of the Fe spin) before a lowest temperature increase 
toward a constant value. For comparison, the same calcula- 
tion has been performed on the kagome lattice which agrees 
with results of Reimers et al. [23| . The values of ^ are labeled 
on the plot. 



a reduction in the magnitude of the susceptibility. 

Since the J2 > magnetic susceptibility appears 
to more closely resemble the experimental data on 
FeCrAs[3] we will focus on that case here. In the limit 
that ^ goes to zero the inverse magnetic susceptibility 
docs not go to the kagome result as it contains a paramag- 
netic contribution due to the uncoupled Fc spins which 
reduces the inverse magnetic susceptibility to zero in- 
stead of reaching a constant, non-zero value as in the 
strictly distorted kagome planes case (see Fig. H]). The 
precise location of this downturn depends on the rela- 
tive magnitudes of the spins on the Cr sites and the Fe 
trimers, here taken to be of unit magnitude. For small 
values of the coupling constant (0 < ^ << 1) between 
the Fe and Cr layers the inverse magnetic susceptibility 
has features reminiscent of a lattice of distorted kagome 
planes, exhibiting a weak upturn followed by a decrease 
at low temperature. For small values of the in-plane cou- 
pling constant (^ >> 1) the inverse magnetic suscep- 
tibility displays properties of both limits. The inverse 
magnetic susceptibility of these intermediate values re- 
duces towards zero as in the cases with small coupling 
within the Cr planes (-^ >> 1), prior to increasing in 



a soft peak at lower temperatures for 



1, before re- 
turning to zero. 

To make contact with the experimental results of Wu 
et al. 18|, it is tempting to choose magnetic couplings Ji 
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FIG. 5: (Color online) A plot of the magnetic susceptibility 
of the Monte Carlo results for -jj- = ^ and -jj- = -j (inset). 



and J2 in such a way as to locate the maximum of the 
susceptibility roughly in the middle of our temperature 
scale as seen experimentally on a scale from to 300K. 
While experimentally FeCrAs is seen to have an easy axis 
below the magnetic ordering temperature, which has not 
been included in our model, we see from Fig. [5] that the 
qualitative features of the magnetic susceptibility seem to 
be reasonably well captured by -j^ ~ ^ when compared 



to Fig. 2 of Ref. [18[- The magnetic susceptibility of 
the Monte Carlo simulation appears qualitatively to be 
an average of the measured magnetic susceptibility along 
the two directions of the FeCrAs lattice. This comparison 
could be made quantitative by converting from dimen- 
sionless units to emu/ (mole Oe) and including magnetic 
moments of larger relative magnitude on the Cr sites than 
the Fe trimers. While for frustrated magnets one often 
finds the inverse magnetic susceptibility to show a Curie- 
Weiss like straight line to extremely low temperatures, it 
is interesting to note that the presence of the Fe trimers 
appears to cause curvature near the ordering transition 
temperature. As such, it may not be necessary to sub- 
tract a temperature independent contribution to obtain 
the straightest possible (x ~ Xo)^^ as carried out by Wu 
et al.^. 



DISCUSSION AND SUMMARY 

It appears that a classical spin model on the layered al- 
ternating distorted kagomc and triangular lattice is able 
to capture magnetic order at a finite temperature with an 
ordering wavevector consistent with that observed exper- 
imentally in FeCrAs. Magnetic susceptibility simulations 
have features qualitatively consistent with measurements 
of Wu et aZ.[18[ for ^ ~ | if we assume unit spins on 
the Cr and Fe trimer sites. It is anticipated that a sim- 
ilar model with quantum spins on this lattice including 



higher orders of magnetic exchange interactions such as 
a Ji, J2, and J3 model where J3 couples the Fe trimers 
may capture the low temperature heat capacity of this 
unusual non-Fermi liquid material. 

It is also plausible that the low temperature phase 
seen in FeCrAs, featuring a Fermi liquid- like constant 
specific heat coefficient coinciding with insulating-like 
(power-law) electronic transport, indicates the presence 
of a U(l) spin liquid insulator with a finite spinon Fermi 
surface as first discussed by Podolsky et at in the context 
of Na4lr308mi- Within such a picture, charge neutral 
spinous lead to gapless excitations and the non- Fermi liq- 
uid behavior of transport may imply that FeCrAs is close 
to a quantum phase transition between metallic and in- 
sulating states. Such a possibility arising in the context 
of the model here presented of stacked distorted kagome 
and triangular planes will be presented elsewhere, along 
with a microscopic derivation of the model [l^. We em- 
phasize that the magnetic ordering obtained in this clas- 
sical Hcisenbcrg model is a robust feature independent of 
whether spinous or electrons are the proper description 
of the elementary magnetic units. Therefore when the 
moment size and/or orbital overlap with the Cr atoms 
decreases, the magnetic transition temperature should 
decrease, as the coupling between Fe and Cr is a crucial 
mechanism for the magnetic ordering. Such physics may 
be accessible in these materials under the application of 
uniaxial pressure. 

In summary, we have studied a classical Heisenberg 
model on the hexagonal Fe2P/ZrAlNi structure. We have 
focussed on the common situation where both the 3f and 
3g sites hold magnetic atoms, with the 3g sites forming 
layers of distorted kagome planes bisected by triangular 
lattice planes. The 3f sites form trimers, small equilateral 
triangles, about each site of the triangular lattice in the 
intermediate planes and in our simple model have been 
treated as single sites of unit spin, as have each of the 
3g sites. The model shows two distinct phase transitions 
depending on the magnitude of the ratio of the coupling 
constants For values of < 1 the 3g spins are 120° 
rotated on the corners of the equilateral triangles, while 
the coplanar 3f spins align opposite(parallel) to the spin 
sum on the nearest 3g distorted hexagon. For > 2 the 
spins take on a ferri(ferro)magnetic form, ferromagnetic 
within each plane. For intermediate values sequential 
transitions are seen with a ferri(ferro)magnetic transi- 
tion at high temperatures preceding a lower temperature 
120° rotation of some component of the canted spin. The 
ordering temperature has been extracted from the heat 
capacity plots for < < inf and a phase diagram 
has been presented, the validity of which in the thermo- 
dynamic limit has been checked by the calculation of a 
Binder cumulant. These results may be of relevance to 
FeCrAs with Cr atoms in the 3g and Fe in the 3f sites. 
We have calculated the magnetic susceptibility for a ra- 
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tio of the coupling constants of ^ = ^ was found to be 
qualitatively similar to the experimental results of Wu et 
al.^. 
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